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We present a Fokker-Planck description of supercooled colloidal systems exhibiting slow relaxation 
dynamics. By assuming the existence of a local quasi-equilibrium state during the relaxation of 
. the system, we derive a non-Markovian Fokker-Planck equation for the non- stationary conditional 

o 



o 



probability. A generalized Stokes-Einstein relation containing the temperature of the system at 



^Nj . local quasi-equilibrium instead of the temperature of the bath is obtained. Our results explain 

Oh. experiments showing that the diffusion coefficient is not proportional to the inverse of the effective 

D . 

jy^ , viscosity at frequencies related to the diffusion time scale. 

I. INTRODUCTION 

> 

The dynamics of slow relaxation systems exhibits peculiar characteristics which make it essentially different from 
the dynamics of systems relaxing in much shorter time scales [J- [(J- The aging behaviour of the correlation functions 
Q and the violation of the fluctuation-dissipation theorem and of the Stokes-Einstein relation [2( are among those 



On 
O . 

significant features which have attracted the interest of many researchers in the last few years. Relaxation phenom- 
ena in glasses, polymers, colloids and granular matter provide innumerable situations demanding new theoretical 

:3 



developments whose implementation constitutes a challenge for nonequilibrium statistical mechanics theories. 



The origin of that peculiar behaviour lies in the fact that during its evolution, the system rests permanently away 
from equilibrium. This feature would justify why results obtained when the system is close to equilibrium are not 
^ \ necessary valid when that condition is not fulfilled. A simple example illustrating this point is supplied from the slow 
, dynamics of a two-level system. During the transition, the system does not thermally equilibrate with the bath and as 
a consequence equipartition law is not valid. It was shown in 10] that coarsening the level of description of the path 
that the system follows in the configuration space going from diffusion to activation regimes leads to violation of the 
fluctuation-dissipation theorem. The fast decay of the diffusion modes makes the system always close to equilibrium 
whereas the activation process occurring at longer time scales does not guarantee validity of the theorem. 

In this paper we will present an example illustrating the peculiar dynamics of slow relaxation systems. Using a 
Fokker-Planck description we derive a generalized Stokes-Einstein relation for supercooled colloidal liquids showing 
that diffusion coefficient and inverse of viscosity are not proportional to each other at low frequencies, when the system 
becomes activated. At higher frequencies, the particles undergo diffusion in the solvent and that relation holds. 

The paper will be distributed as follows. In Section II, we will briefly review the generalization of the Onsager 
fluctuation theory to nonequilibrium aging states. Section III will be devoted to the formulation of a mean field 
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theory describing the dynamics of high concentrated suspensions. This theory will be used in Section IV to interpret 
experimental results supporting the violation of the Stokes-Einstein relation in supercooled colloidal liquids. In the 
conclusions section we will summarize our main results. 

II. SLOW RELAXATION DYNAMICS AND LOCAL QUASI-EQUILIBRIUM 

When the relaxation of a system is slow enough, its evolution takes place through a sequence of states characterized 
by a set of thermodynamical quantities in which the intensive parameters are, in general, different from those of the 
bath. It is then said that the system is at quasi-equilibrium. 

A statistical mechanics description of those systems involves the non-stationary conditional probability density 
P(a ,to\a, t), with a — (ai, ■■■,ot n ) representing the state vector of the system, and a the initial state at time to. 
The local quasi-equilibrium states are characterized by the probability density P qe (a,t) satisfying 



lim P(a 0) t |ai,i) = Pqe{a,t). (1) 

t — >oo 

When relaxation takes place in a short time scale, P qe (a, t) reduces to the local equilibrium probability density in 
which the intensive parameters are those of the bath. 

The evolution of the conditional probability density is governed by a Fokker-Planck equation. To obtain it, we will 
first formulate the continuity equation 

d d \ 1 

— P(a ,to\a,t) = -— ■ |P(So>*o|a,%«(a Q ,a;to,£)J , (2) 

in which Pv a represents the phase space probability current with v a the stream velocity in a-space. That velocity 
can be obtained from a mesoscopic version of the noncquilibrium thermodynamics for slow relaxation systems we have 
proposed [Tl| . 

The existence of a local quasi-equilibrium state in a-space enables one to formulate the Gibbs equation 

T(t)5s qe (t) = Se qe (t) - J v qe {t)SP{t)da, (3) 

where s qe (t) is the entropy and e qe (t) the mean internal energy per unit mass at quasi-equilibrium . We will assume 
that the temperature T(t) of the system at local quasi-equilibrium is only a function of time whereas the chemical 
potential n qe (t) can be a function of a. 

The entropy of the system can be expressed in terms of the probability density by means of the Gibbs entropy 
postulate [13 in the form 



s(t) = ~J P( to \t)ln^lda + Sqe (t), (4) 

where k is the Boltzmann constant and m the molecular mass. 

The form of the current Pv a can be inferred from the entropy production by relating it, as done in noncquilibrium 
thermodynamics |l2j . with the corresponding thermodynamic force. The entropy production is obtained by combining 
the rate of change of s qe (t) obtained from Eq. © with the time derivative of Eq. QJ. After using Eq. and 
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integrating by parts assuming that P (iol^)^ct(^Oj t) vanishes at the boundary, one may identify the non-equilibrium 
chemical potential /i(to;t) = fJ- qe {t) + fc ^^ In P p to ^ and the corresponding thermodynamic force , which is 
conjugated to the probability current Pv a . Following the scheme of nonequilibrium thermodynamics one simply 
obtains 



nJto\t) = -R{t ;t) 

where we have defined the generalized force 



X(a,t) + ^-?-lnP(t \t) 
m oa 



(5) 



X(a,t) = - — 
oa 



kT{t) 



\nP qe {t) 



(6) 



The time dependent transport coefficients B_(to ; t) (related to the Onsager coefficients) incorporate memory effects 
through its dependence on time jl^j . and the local quasi-equilibrium probability density is given by 



P qe (t) = p e e*9m^°- e «*\ (7) 

where P e is the equilibrium distribution function. Eq. (J7J) can be obtained by using the statistical definition of the 
entropy in the Gibbs equation © , Q| • Now, by substituting Eq. into (J2J the resulting Fokker-Planck equation is 



dP(t \t) d 



Of 



dc 



B(t ;t) 



X(a, t)P(t \t) 



kT(t) dP(t \t) 
m da 



(8) 



Using the Fokker-Planck equation JSJ , we can calculate the evolution equation for the equal-time correlation function 
of a- variables which for sufficiently long times leads to the following expression for the temperature of the system 



T{t) = j(X(t)-a(t)) qe , (9) 

with (...) de noting average at local quasi-equilibrium. This result reduces to the one corresponding to fast relaxation 
processes Wa in which the average is performed at local equilibrium. 

In Ref. [llj, we have shown that the temperatures of the system and of the bath Ts{t) are related by T(t) = ATsit), 
expressing lack of thermal equilibration. Equilibrium is reached when A = 1 in which case, for linear thermodynamic 
forces, one recovers the expression for equipartition law. The presence of T(t) in the last term of Eq. Q, gives rise 
to a modified version of the fluctuation-dissipation theorem. 



III. FOKKER-PLANCK DESCRIPTION OF CONCENTRATED COLLOIDAL SUSPENSIONS 

In this section, we will study diffusion in concentrated colloidal suspensions by means of an effective medium 
approach elaborated on the grounds of the theory we have introduced in the previous section. We will analyze the 
motion of a test colloidal particle through the suspension by using the Fokker-Planck equation describing the evolution 
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of the two-time probability density P(ctQ, to\a, t). In this case, the phase space vector is simply a = (u, r), where u 
and f represent the velocity and position of the test particle. 

To derive the Fokker-Planck equation, we will first formulate the Gibbs equation 

TpSs = p5e — p~ l p8p — ml pSPdu, (10) 



which incorporates the effects of interactions among particles through the term containing the excess of osmotic 
pressure p. Here e is the internal energy, m the mass of the particle and p — m J Pdu the mass density. In this case, 
the quasi-equilibrium probability density is given by 

Here, the interaction of the test particle with the other colloidal particles is represented by means of the term p p 
which can be expressed in terms of a virial expansion in p, whereas p qe constitutes the ideal chemical potential. 
Introducing the fugacity z — Pa, with the activity coefficient given by a = e^r p p , Eq. can be expressed as 

kn f z , . 

s = - Pin — du + s qe , (12) 

P J Zqe 

where z qe — e^^ qa ~^ u 1 is the fugacity at quasi-equilibrium. Notice that in the limit of infinite dilution p~ 1 p = 0, 
thus implying a = 1, which corresponds to the ideal case. 

The existence of interactions among the test and the other colloidal particles, gives rise to a force term in the 
continuity equation for P. In particular, according to Eq. 1)11(1. the force entering into this term is given by — V(p~ 1 p) 
(see for example Ref. 17]). Now, by taking into account the definition of the activity coefficient a, the continuity 
equation for P becomes 

~kT(t) 



dP 

— + V • Pu - V 
at 



In a 



dP d 

^ = -^-Pvu, (13) 
ou ou 



rn 

where Pv^ is the probability current, with the stream velocity in u-space. 

In accordance with the general formalism of nonequilibrium thermodynamics, the expression of the probability 
current Pv^ can be obtained from the entropy production which only contains dissipative terms and is simply given 
by 

— <»» 

Eq. i|14fl has been obtained by taking the time derivative of Eq. (|12|) . using Eqs. (|13fl . (|llfl and JTUJl, and the balance 
equation for the energy e(t), neglecting viscous dissipation [12|. In order to obtain Eq. (|14|l . we have also identified 
the nonequilibrium chemical potential 

p(t ;t) = p qe (t) + ^\n^l. (15) 
m z qe it) 
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FIG. 1: Comparison between the viscosity and diffusion coefficients for two vaiues of the voiume fraction. Triangles and squares 
represent viscosity data from Ref. The gap between viscosity and diffusion at short frequencies is proportional to the scaling 
factor A of the temperature. At larger frequencies A — * 1. 

By taking into account the definition ©, the resulting generalized Fokker-Planck equation for the two-time prob- 
ability density is 



dP „ 

^7 + V ' UP 
at 



kT(t) 



In a 



cm au 



(16) 



where f3 characterizes the dissipation of the kinetic energy of the test particle. This coefficient is, in general, function 
of time, position and density; its dependence on time is a consequence of the non-Markovian nature of the stochastic 
process Q, Q, Q. 

FromEq. JT§J}, we will now derive the evolution equation for the velocity field by taking the time derivative of the 
momentum pv = m J uPdu and using Eq. (|16|) . After integrating by parts one obtains 



dv n _ kT(t) 

p— = -Ppv 

at m 



d In a 
<91np 



(17) 



where we have used the fact that, in the absence of an externally imposed How, the second moment P = m J (%t — 
v)(u — v)Pdu can be approximated by P = — pi, |is| . 

At times such that t 3> the particle enters the diffusion regime in which one may neglect the time derivative 
in Eq. H17|) . obtaining then the constitutive relation for the mass diffusion current pv which substituted into the mass 



continuity equation 



di 



-V • (pv), yields the Smoluchowski equation 



(18) 
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From Eq. (|17fl we may identify the diffusion coefficient 



D 



kT(t) 

m(3 



1 



d In a 
dhip 



(19) 



The form of Eq. I|19|) coincides with the expression given in Ref. |l2|. Notice, however, that it contains the temperature 
of the system at quasi-equilibrium instead of the temperature of the bath. 



A. The generalized Stokes-Einstein relation in supercooled colloidal systems 



It has been shown in |9j that in supercooled colloidal liquids at short frequencies the Stokes-Einstein relation is not 
fulfilled. Measurements of the viscosity and the diffusion coefficients obtained in Refs. [2^, indicate the existence 
of a frequency domain in which those quantities are not proportional (see figure I). It is our purpose to offer here an 
explanation for those experimental results. 

In [3| it was shown that the colloidal system can be assumed isotropic and homogeneous at the considered volume 
fractions, ranging from the freezing to the glass transition values <f> g . Consequently, the transport coefficients can be 
assumed to be independent of the wave vector k. We will first introduce the effective viscosity 

kT 

= ^D-\ (20) 

where R is the radius of the particle. 

The dependence of the viscosity on the frequency can be inferred through the expression of the penetration length 
A of the diffusive modes. In the dilute case, A is given through A -1 = (™) 5 , with r = -^j- the characteristic 
diffusion time [22| and Do = g~^jb the diffusivity in the pure solvent, with 770 its viscosity. At higher concentrations, 
it can be assumed to be 



A- i = f ^(™r, (2i) 

where the exponent 7 and the scaling function T may, in general, depend on the volume fraction. The quantity TR 
can be interpreted as an effective radius. 

The penetration length gives an estimation for the average size of a cage formed by the surrounding particles 
At high frequencies, the test particle performs a free Brownian motion inside the cage and the viscosity is that of the 
pure solvent. At short frequencies, collisions of the test particle with the other particles modify the viscosity. Thus, 
we will assume the following expression for the effective viscosity 

ri((j) = vo (1 + ^A- 1 ) , (22) 

valid up to first order in RX^ 1 reflecting the fact that cage diffusion is important when R ~ A, [2^ . 

By taking into account the general relation between system and bath temperatures obtained from our formalism, 
T = AT B , PJ, Eq. igDJl can be written as 
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which constitutes a generalization of the Stokes-Einstein relation to systems exhibiting slow relaxation dynamics. 

From Eq. and from its definition, the quantity A is essentially proportional to the velocity moments of the 
particle. In particular, when X(t) is a linear function, A is proportional to the velocity correlation function. It is 
then plausible to assume the following expansion 13 



A= [l+6(rw) e ] 



(24) 



where b and the exponent e may, in general, be functions of the volume fraction <f>. Comparison of Eq. (|24J) with 
the corresponding relation obtained in Q by fitting the experimental data yields In b = e (40^ - 37^ + 2.659, and 
e = -0.77. 

To contrast the theoretical expression of the effective viscosity (1221 with the experimental results reported in 
we will take the logarithm of Eq. (|22|l and expand the result in terms of the variable logrw, around tuj — 1. Up to 
first order we obtain 



log 77 = log 



rio 1 + 



i + r 



logrw. 



(25) 



A direct comparison of Eq. (|25|l with the fitting of experimental data given by the expression log 77 = cq + s(<ft) log rw, 
with s(<p) « 6.39 — 13<^> and cq a constant, gives the dependence of 7 and T on the volume fraction. 



8 



The expression of the diffusion coefficient at the frequencies we are considering can now be derived by using Eq. 
(|23|l . In figure 2, we represent the diffusion coefficient given through Eq. (|23[1 as a function of the reduced frequency. 
At higher concentrations (<p — 0.56, 0.55), a closer agreement with the experimental values of D given in [jj has been 
obtained with e = —0.77 (solid lines). At lower concentrations (</> = 0.53, 0.52), a closer agreement with experiments 
has been obtained with e = —0.67 (dotted lines). This dependence of the fitting on the value of the exponent e 
suggests that it is also a function of the volume fraction. 

IV. CONCLUSIONS 

In this paper, we have proposed a general formalism to analyze the dynamics of systems which relax in long time 
scales of the order of the observation time or even longer. The fact that the system always evolves through transient 
states never reaching equilibration with the bath is the origin of a peculiar behaviour different from that of systems 
whose propagating modes relax rapidly ending up in a quiescent state. The dynamics is non-Markovian and the 
correlation functions exhibit aging effects. 

We have formulated a generalization of Onsager's theory to nonequilibrium aging states to derive a Fokker-Planck 
equation which captures the main characteristics of the dynamics such as non-stationarity through the two-time 
probability density, non-Markovianity through the time-dependence of the transport coefficients and lack of thermal 
equilibration between system and bath. This equation is the main result of our analysis from which the behaviour of 
the correlation functions follows. 

We have applied the theory to study the relaxation in supercooled colloidal liquids for which violation of the Stokes- 
Einstein relation at low frequencies has recently been found experimentally. The relation between diffusion coefficient 
and viscosity we have obtained explains those experiments. The key ingredients in this interpretation are the existence 
of a quasi-equilibrium state with a temperature different from that of the bath which is proportional to the velocity 
moments and the fact that the viscosity depends on the penetration length of the diffusion modes resulting from the 
underlying activated process taking place at sufficiently long times. 

The formalism proposed thus offers an interesting theoretical framework for the study of glassy behaviour in colloidal 
liquids. 
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